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Abstract. It is conjectured that every integer N > 454 is the sum 
of seven nonnegative cubes. We prove the conjecture when iV is a 
muhiple of 4. 

1 Introduction 

Waring famously asserted in his Meditationes Algebraicae of 17700 

Omnis integer numerus vel est cuhus, vel e duobus, tribus, 4j 5, 6, 1, 8, vel 
novem cubis compositus . . . 

meaning that every positive integer is the sum of at most nine positive cubes 
(equivalently, of exactly nine nonnegative cubes). A proof of this was given by 
Wieferich in 1909 (with an error later patched by Kempner [1912]). Landau 
then proved that only a finite set requires nine cubes, and Dickson [1939] iden- 
tified this set with {23,239}. Linnik then established, ineffectively, that only a 
finite set requires eight cubes, and seven suffice after some point [Linnik 1943] . 
In 1984, McCurley gave an effective proof of Linnik's result, demonstrating 
that every integer larger than exp(exp(13.94)) is the sum of seven positive 
cubes [McCurley 1984| . This was recently reduced to exp(524) [Ramare 2007) 
using an analytic sieve argument. 

It is believed that the exceptional set for Linnik's seven cubes theorem is 

{15, 22, 23, 50, 114, 167, 175, 186, 212, (1) 
231, 238, 239, 303, 364, 420, 428, 454}. 

(Indeed one expects that every sufficiently large integer is the sum of four pos- 
itive cubes |Deshouillers et al. 2000| . but even such a statement with four re- 
placed by six is well beyond our ability to prove.) 

In the other direction, it is shown in [Bert ault-Ramar e-Z immer mann 19991 The- 
orem 1 and Lemma 3] that if 454 < N < 2.5 • 10^^ then N is the sum of cubes 
of seven nonnegative integers; and it is observed in [Ramare" "20071 p.60] that 



^ Ellison [1971, p. 10] reports that this statement appears on pages 203-204 of the 1770 
edition. In the EngUsh translation Waring 1782 of the 1782 edition, this statement appears 
on page 336 as part (9) of "Theorem 9"; see also the discussion on page 379 of the same 
translation, and the first section of Ram are 20071 for a fuller treatment of the history of this 
problem than we give here. 
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the computation reported in |Deshouillers et al. 20001 433-434] raises the up- 
per bound to exp(78.7) > 10^"'. But the computation to raise this bound to 
exp(524), and thus prove that H]) is the full exceptional set, remains utterly 
infeasible. 

We give a different kind of partial result, where N is restricted by a congruence 
condition but not by size: 

Theorem. Every multiple of 4 except 212, 364, 420, and 428 can be written as 
the sum of seven nonnegative cubes. 

This is not the first such result, but the only earlier work in this direction that we 
know of is the proof in |Bertault-Ramare-Zimmermann 1999j that if = or 
±1 mod 9 and N is an invertible cubic residue mod 37 then N is the sum of seven 
nonnegative cubes. Bertault et al. note that 37 could be replaced by various 
larger primes congruent to 1 mod 3. But the condition mod 9 is essential, and 
restricts mod 9 to the three most common residues for a sum of seven cubesH 
Compared with [Bertault-Ramare-Zimmermann 1999] . our new ingredient is 
the use of a quadratic form Q = X]i=i '^i^f ^it^ (01,02,03) 7^ (1, 1, 1) which 
is nevertheless known to represent all positive integers in certain arithmetic 
progressions. 

The rest of the paper is organized as follows. We review the basic identity ([4]) 
for writing suitable integers as sums of six integer cubes, then give a crite- 
rion (Proposition 1) under which the cubes are positive. The criterion re- 
quires an auxiliary prime p = 2 mod 3 in an interval (AN-^^^ , BN^^^) and a 
small enough positive integer xq such that p\N — Xq. We then choose Ci, C2, C3 
and show that a suitable Xq exists provided p satisfies a congruence condition 
mod 72. Finally we use the explicit bounds of [Ramare-Rumely 1996| on the 
distribution of primes in arithmetic progressions, plus a short further compu- 
tation of prime chains, to prove that such p exists if > A'^o 10^^. This 
completes the proof because that Nq is well below the threshold of 2.5 • 10^^ 
of [Bertault-Ramare-Zimmermann 1999j . 

Since we formulate the proof so as to use only positive rather than nonnegative 
cubes, we automatically get a representation of A^ as a sum of seven positive 
cubes for all A^ > 10^^ divisible by 4. For A^ < 10^^, there are cases not listed 
in ([1]) for which A^ is a sum of six or fewer nonnegative cubes but not exactly 
seven; the largest of these is apparently 2408 (as it happens a multiple of 4), 

^ The following table gives the distribution mod 9 oi N = + n'^ + . . . + among the 
9^ possibilities of {rii , 712 , • ■ ■ , '^7) mod 9: 



N mod 9 





±1 


±2 


±3 


±4 


proportion 


17.97% 


16.32% 


12.21% 


7.68% 


4.80% 



Note that of the seventeen exceptions listed in JTJ, eleven are congruent to ±4 mod 9, and the 
remaining six to —3 mod 9. These rare residues ±3, ±4 are also the least easily accessible to 
current approaches to the seven-cube problem, including ours: it will be seen that we must 
work hardest to prove our theorem for N congruent to ±3 or ±4 mod 9. Note that while our 
condition 4\N also puts in a probabilistically favored congruence class, the discrepancy is 
minuscule: a random sum of seven cubes is divisible by 4 with probability only 25.39%. 
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whose only representations as a sum of seven or fewer positive cubes are 



2408 10'^ + 10^ + 7^ + 4^ + 1^ 

= 12^ + 8^ + 5^ + 3^ + 23 + 23 

= 113 + 83 + 63 + 63 + 53 + 23 

= 103 + 103 + 63 + 43 + 43 + 43 (2) 

and permutations of these four sums. 
2 A six cube identity 

The use of the identities equivalent to the following Lemma and its Corollary 
to study sums of seven cubes goes back at least as far as jLinnik 19431 §2]. 

Lemma 1. Let 

Q = Ec.Xf, C = ^c3. (3) 

1=1 1=1 

Then 

3 

{{c^P + + (c,p - X,f) = 2Cp3 + QpQ. (4) 



Proof: Apply the identity [r + s)3 + (r — s)3 = 2r3 + Qrs^ to each of the three 
terms in the sum. □ 

Corollary. // p and the Ci are positive integers and the Xi are integers such 
that \Xi\ < Cip for each i, then 2Cp^ + &pQ is the sum of cubes of six positive 
integers. □ 



3 Strategy 

Given N , we choose integers ci > and a prime p = 2 mod 3 of size roughly 
N^/^, and let a;o be a small nonnegative integer such that N — Xq — 2Cp'^ + 6pQo 
for some integer Qq > 0. Such an xq exists if N/p'^ is large enough because 
every integer is a cube mod 6p. If, moreover, xq can be chosen such that Qo is 
represented by the quadratic form Q, and N/p^ is small enough so that all the 
terms in the sum in the identity ^ are positive, then we can use that identity 
to write TV as the sum of cubes of seven positive integers. 

Because no ternary quadratic form represents all Qo > 0, we may need to put 
a further congruence condition on xq modulo some P relatively prime to 6p; we 
shall then choose the least xq > satisfying these congruences, so that xq < 6f3p. 
We find that this imposes the following lower and upper bounds on N/p^: 

Proposition 1 Let ci, 02,03 be positive integers with c\ = minci. Set C = 
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Y^'l=i • ^or some (3 > 1 assume that 

2C + 2I6/33 < 4 < 2C + 6c?. (5) 
pi 

If xq is a positive integer such that xq < 6/3p and 

N - xl = 2Cp^ + 6pQo, (6) 

then Qo > 0; if further 

Qo = c,Xf + c^Xi + c^Xl (7) 

for some integers Xi, then N is the sum of cubes of seven positive integers. 
Proof: The lower bound on N /p^ assures that 

N -xl>N - {<ol3pf > 2Cp^ 

so Qo > 0. Given a solution of ([7]), we use the identity ^ to write N — Xq as 
the sum of cubes of six integers. Since xq > 0, it thus suffices to verify that 
< dp. Indeed we have 

2 N^xl-2Cp^ N - 2Cp^ 6c?p3 32.32 

QAj <Qo^ < < — — = < c,p 

bp bp bp 

so Xf < cfp^ = [ciPi)"^. Since Cip > 0, we are done. □ 
The inequalities Q require 

/3 < ci/s/se. (8) 

Assuming this condition holds, ^ restricts p to an interval {AN^^^^, BN^^"^) for 
some constants A, B with < A < B. We then use explicit bounds for the 
distribution of primes in arithmetic progressions to find A'o < 10^® such that a 
suitable p exists for all N > Nq. 

When is a multiple of 4 but not of 8, we shall need to impose a condition on 
p mod 8; when N falls in one of the hardest residue classes ±4 mod 9, we shall 
also impose a condition on p mod 9. 



4 Choices and analysis 

We choose (ci,C2,C3) = (4/3,4/3,6/3) where /? = 1 or /? = 5 depending on the 
residue of N mod 9 (as specified later in this section). Then condition ([8]) is 
satisfied, and 

CiX2 + c,X^ + c,Xl ^ 2fi{2Xl + 2X1 + 3^1) • (9) 

We calculate C = 344/3^, whence A = /3- 7^^1072 and B = /3^Vv^, with 
ratio BjA = (134/113)^/^ > 1.0584. 
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We choose xo so that 

xl = N- 2C/ = iV - 688(/3p)^ mod 6^p, (10) 

which is possible because every integer is a cube mod 6(3p. We select the least 
positive xq satisfying this congruence; thus xq < 6/3p. As N and 6(3p are even, 
so is Xq. Since also 4|7V and 2|C, while (3 and p are odd, it follows that N—2Cp^ 
is a multiple of 12/3p. That is, 

TV - x§ - 688(/3p)-^ 
12^^ 

is an integer. 

This integer Qi is positive by our Proposition. Set Qo = 2(3Qi. In view of 
we need Qi to be represented by the quadratic form 2Xi + 2X| + 3X|. This 
quadratic form is unique in its genus, so it represents all nonnegative integers 
that are not excluded by congruence conditions. In this case this means all 
Qi that are neither congruent to 1 mod 8 nor of the form 9* (9m + 6) for some 
nonnegative integers m,t. (See [Dickson 1927| (16), pages 44-45] for this char- 
acterization of the integers represented by Qi.) It remains to choose p so that 
these conditions are satisfied. 

The condition for Qi mod 8 holds automatically if = mod 8, because then 
the numerator N — xf^ — 688{f3p)^ in (fTT|) is a multiple of 8 while the denominator 
is not, so the quotient Qi is even. Assume then that A^ = 4 mod 8. Since 
688 = 16 mod 32 and /3p is odd, we have 688(/3p)'^ = 16 mod 32 as well. Since 
Xq is even, Xq is a multiple of 8 not congruent to 16 mod 32. Therefore 

4 4 

We therefore choose p so that 

A 

3/3p EE — mod 8, (13) 

and this guarantees that Qi ^ 1 mod 8, so Qi passes the mod-8 test for repre- 
sentability by the quadratic form 2A^ + 2A| + 3A|. 

Next we ensure that Qi ^ (9m + 6)9* by choosing (3 G {1, 5} and p mod 9 so 
that Qi ^ 0, 6 mod 9. Since p = 2 mod 3 implies p^ = —1 mod 9, the choice of 
P determines 688(/3p)^ mod 9. Thus, as 688 = 4 mod 9, we have 

A^ - 688(/3p)3 = A + 4/3^ mod 9. (14) 

Now Xq = or ±1 mod 9 for all integers xq. Therefore if 

A + 4/5=^ ^ or ±1 mod 9 (15) 

then N — Xq — Q88{/3p)^ cannot be a multiple of 9, whence Qi is not a multiple 
of 3, so a fortiori not congruent to or 6 mod 9. We have = 1 mod 9 for 
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f3 — 1, and = —1 mod 9 for /3 = 5; thus, unless N = ±4 mod 9, we may 
choose (3 so that N satisfies condition ([15]) 

In the remaining cases N = ±4 mod 9, we choose /3 so that 3\xo by requiring 
^ = 5 for AT EE 4 mod 9 and /? = 1 for iV = -4 mod 9. Then 27\xl, and the 
numerator N — Xq — 688{(3p)^ in pT]) is a multiple of 9 that we can control mod 27 
by choosing p mod 9. Indeed if /3p = 3fc ± 1 then {(3p)^ = 9fc ± 1 mod 27. Since 
688 = 13 mod 27 this gives 

N~xl- 688(/3p)3 TV ±13 

^ = g k mod 3, (16) 

with the sign chosen so that 9\N ± 13. Dividing by 4/3p, we conclude that 

^l^p(k-^:i±l^)mod3. (17) 



3 V 9 , 

We thus choose A: ee /3 + ((iV ± 13)/9) mod 3. That is, 



/3p = 3/3 + ^"t^^ ± 1 mod 9, (18) 



with the sign in ±1 chosen so that ±1 = —(3 mod 3. Then Qi/3 = 1 mod 3, so 

^3- 



Qi passes the mod-9 test for representability by 2Xf + 2X| + 3X| 



5 Conclusion 

To finish the proof of our Theorem, we show: 

Lemma 2. For (3 £ {1,5} let A = p-^/y^WH and B = / ^/ml. Set 
No ~ 10^^. Then whenever N > Nq there exists a prime p E 
in each odd congruence class I mod 72 = 8 • 9 with I = 2 mod 3. 

This will suffice because 10^^ is smaller than the lower bound of 2.5 • 10^^ 
of [Bertault-Ramare-Zimmermann 1999] on an integer iV not listed in ([1]) that 
is not the sum of seven nonnegative cubes. 

Proof: Taking fc = 72 in [R amare-Rumely 1996] Theorem 1], we find that for 
every prime p > 10^^ there exists a prime p' > p such that p' = p mod 72 
and p' < ((1 + £72)7(1 — £72))^- Consulting [Ramare-Rumely 1996 §5, p. 419, 



Table 1], we find £72 < 0.013, so (1 + £72)7(1 — £72) < 1.027 which is well below 
the gap ratio of B/A = (1347113)^/^ > 1.0584 that we need. This estabhshes 
Lemma 2 for TV > 1072(5 • 10^°)^. 

That bound is not small enough for our application because it exceeds the 
threshold of [Bert ault-Ramar e-Z immer mann 1 9 9 9 j and even (by a factor of 13.4) 
the improved bound of lO^'' reported in jRamare 20071 p. 60]. But it reduces the 
proof of Lemma 2 to a finite computation. To make this computation manage- 
able, we prove: 

Sublemma. In each congruence class I mod 72 coprime to 72 there exist an 
integer Mi and primes Pi (0 < i < Mi) such that po < 19541, pi-i < pi < 
1.0584pj_i for each i = 1, 2, . . . , Mi, and pj^j^ > 10^". 
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To derive Lemma 2 from this Sublemma, let p be the smallest prime such that 
p = I mod 72 and p > AN^/^. Then p > 19541 because 

N > No = 10^^ > 1072(5 • 19541)^ 

If p > 10^", use [Ramare-Rumely 1996| . Else apply the Sublemma and find the 
maximal nonnegative i < Mi such that p > pi. Then 

p < K+i < 1.0584 p, < 1.0584 p < 1.0584 AN^^^ < BN^'^, 

and we are done. 

(Conversely, the existence of the gap ratio 17573/16493 < 1.0655 between the 
primes 16493 and 17573 congruent to 5 mod 72 shows that A'o cannot be brought 
much below 10^^ in Lemma 1.) 

To find the primes pi required by the Sublemma, wc run the following algorithm 
for each I with some choice of positive 5 < 0.0584: 

• Let Pa be the largest prime p < 19541 such that p = I mod 72. Set i — 0. 

• While Pi < 10^", let pi+i be the least prime such that pi+i = pi mod 72 
and Pi+i > (1 + 5)pi, and increment i to i + I. 

Each "largest prime" and "least prime" is found by simply stepping down or 
up the congruence class until a prime is found. The second step is repeated 
at most log(10^°)/ log(l + S) ^ 23/6 times. Once pi exceeds 10^*^, we suc- 
ceed if maXj^^j pi^i/pi < 1.0584; otherwise we try again for a smaller S: de- 
creasing 6 lengthens the computation but skips fewer primes. We find that 
6 = 0.01 is small enough for the computation to succeed for each I. Then 
Ml < 1250 in each case, and the largest pi+i/pi ratio occurs at = (5,1), 
namely 21101/19949 < 1.0578 < 1.0584. This computation, programmed in 
PARI-GP [Batut et al. 1998] . takes less than a minute to run on an office desk- 
top machine (G5 PowerMac with a 1.8 GHz processor), and completes the proof 
of the Sublemma and thus of our Theorem. 

Remark: We could also have used [Ramare-Rumely 1996[ Corollary 5.2.2] to 
prove a result similar to Lemma 2 but with the weaker bound iVo = 10^"'^, 
which is still good enough because 10^^ < 2.5 • 10^®. But this would imphcitly 
rely on the much more extensive calculations by Ramare and Rumely of the 
distribution in arithmetic progressions of the primes up to lO^'', which involve 
the computation of hundreds of millions of primes, whereas we needed fewer 
than 1250 (/i(72) = 3 • lO"' primes to prove Sublemma 2. 
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